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Huygens' principle following from the d'Alembert wave equation is not valid in two-dimensional 
space. A Schrodinger particle of vanishing angular momentum moving freely in two dimensions 
experiences an attractive force - the quantum anti-centrifugal force - towards its center. We connect 
these two phenomena by comparing and contrasting the radial propagators of the d'Alembert wave 

■ equation and of a free non-relativistic quantum mechanical particle in two and three dimensions. 

I. INTRODUCTION 

< 

In the year 1917 Paul Ehrenfest |1[ addressed the question: Why is the space we live in three-dimensional? He 
provided three physical phenomena known at his time that crucially depend on the number of space dimensions: (i) 
The generalized electrostatic attraction leads only in two and three dimensions to stable circular orbits, (ii) only in 
three dimensions do we find as many electric as magnetic field components and (iii) Huygens' principle, based on the 
d'Alembert wave equation, is only valid in a world with an odd number of spatial dimensions. 
I Iii the present paper we connect this phenomenon of the breakdown of Huygens' principle in two dimensions to 
the surprising appearance of an attractive force acting on a free particle of vanishing angular momentum in two 
dimensions Q2J. We show that both phenomena have the same physical origin: Interference of waves. 

We compare and contrast the Green's functions of the free Schrodinger equation corresponding to vanishing angular 
momentum in two and three dimensions. Since the rotational symmetry of the initial wave function is preserved in 
the time evolution, we can confine ourselves to the Green's function <?' 2 ' and g^ of the radial motion. We derive 

■ analytical expressions for these quantities by taking the product of the Green's functions of the individual one- 
i-G | dimensional motion along each Cartesian axis and integrating over the angular part. 

We show that the radial Green's function in three dimensions is the difference of two Green's functions of two free 
■y ' one-dimensional motions: One Green's function describes the free motion from the initial radius p to the radius r. 
The other Green's function corresponds to the motion of the particle starting at the mirror image — p of the initial 
radius p reaching r. In the three-dimensional space the phase difference between these two Green's functions is tt. 
q-i In contrast, the corresponding Green's function of the two-dimensional radial motion is more complicated. It 
is again the difference of two one-dimensional Green's functions. However, they do not correspond to free motion 
anymore. Moreover, the phase difference between the two contributions is now n/2. This feature reflects the fact that 
in two dimensions the radial Schrodinger equation for vanishing angular momentum enjoys an additional potential 
that is attractive and of quantum origin. It corresponds to the quantum anti-centrifugal force. 

Whereas in one and three dimensions this potential is absent, it reappears in four and higher dimensions. However, 
it is repulsive. In contrast to the breakdown of Huygens' principle, which takes place for all even dimensions, this 
additional potential displays the unique behaviour of being attractive only in two dimensions. 

Our article is organized as follows: In section 2 we briefly review the physical origin of the breakdown of Huygens' 
principle and summarize the essential features of the quantum anti-centrifugal potential. Here we confine our dis- 
cussion to two-and three-dimensional situations. In section 3 we derive analytical expressions for the radial Green's 
function of the free Schrodinger equation in two and three dimensions. In the appendix we make contact with the 
corresponding expressions for the d'Alembert wave equation. 



II. TWO DIMENSIONS ARE VERY DIFFERENT FROM THREE 



In the present section we review two phenomena which appear to be unrelated at first sight- the breakdown of 
Huygens' principle applied to d'Alembert waves in two dimensions and the quantum anti-centrifugal force acting on 
a Schrodinger particle moving freely in two dimensions. Here we do not dwell on mathematical formalism but focus 
on the key ideas. For the sake of completeness we include detailed derivations of the propagators of the d'Alembert 
wave equation in the appendix. 
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A. Breakdown of Huygens' principle 



In El Ehrenfest summarized earlier work by J Hadamard on the propagation of a signal according to the d'Alembert 
wave equation. When we express his findings in today's language, he investigated the Green's function of the 
d'Alembert equation: A delta function perturbation orginating at time t = from the origin travels in three di- 
mensions as a spherical wave which at time t had reached r. In this case the perturbation always keeps its shape, 
since for t > the propagator in three dimensions reads 

G rct (*, r) = ^ 47]T ^ 

However, in two dimensions the delta function perturbation does not propagate as a well localized circular wave, but 
contains a tail, a so-called wake field . As shown in the appendix, the propagator takes the form 

G 2(M = e[(c()2 - 21 



^(ci) 2 - r 2 

Due to the Heaviside step function there exists a sharp wave front. Moreover, together with the square root in the 
denominator it creates a long tail of excitation in the inside domain of the propagating circular ring []. This wake 
field can be observed with water waves after a stone has excited a wave in an initially quiet pond. 

It is interesting to trace the mathematical reason for the occurrence of the square root in the propagator back 
to the dimensionality of space. Indeed, we show in the appendix that the square root stems from the fact that the 
two-dimensional volume clement reads drrdip. In contrast the three dimensional volume clement drr 2 d9 sin 6 dip has 
in the angle variable an additional factor sin 8. It is because of this extra factor that the interference of plane waves 
in three dimensions gives rise to a spherical Bessel function 

I r d 9 S in9e-^- e ^ S ^-^j (kr) (1) 



2 Jo kr 
whereas in two dimensions we find an ordinary Bessel function 

7T 

d(pe -ikroos V = J ( fcr ) ( 2 ) 



1 

2^ 



The spherical Bessel function 



Jo(fcr) = — (e ikr - e~* kr ) (3) 



originating from the interference in three dimensions is an exact superposition of an incoming and an outgoing spherical 
wave with a fixed phase difference of n. 

In contrast, the familiar asymptotic expansion Q of the ordinary Bessel function 



Jo(fcr) - v i cos {kr = V ^ e-W4 {elkr - ie ^ kr) (4) 

suggests that the interference of plane waves in two dimensions gives rise to an approximate superposition of an 
incoming and an outgoing cylindrical wave with an asymptotic phase difference of tt/2. However, equation (Jij) is 
only the most elementary asymptotic approximation of the Bessel function. In the literature there exists a variety 
of techniques || to find more accurate expansions. Independent of their specific form, they clearly indicate that the 
result of the interference of plane waves in two dimensions is not a pure circular wave. This also manifests itself in 
the fact that the zeros of Jq are not equally distributed like jo but rather bunch towards the origin H . 



1 In his novel Flatland E Abbott [|j describes the phenomenon that people living in a two-dimensional world have a dim vision 
because of fog. They see unsharp surroundings. We do not know if the author was aware of the mathematical and physical 
details of Green's functions in two dimensions. 
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B. Quantum anti— centrifugal force 



We now turn to a seemingly different topic: The Schrodinger equations of a free particle of mass M and energy 
E > in two and three dimensions. We recall that in two dimensions the radial wave equation for vanishing angular 
momentum (m = 0) contains an attractive potential. This potential is related to interference of plane waves in two 
dimensions as expressed by the Sommerfeld representation (||) of the zeroth Bessel function. In this way we make the 
connection to Huygens' principle. 

We start from the time independent Schrodinger equation 

h 2 

of the free particle and introduce polar coordinates r and ip by the transformation x = r cos tp and y = r sin (p. With 
the help of the ansatz 



V(x,y) = e 



we arrive at the radial Schrodinger equation 

f d 2 2M 

Here we have introduced the effective potential (t]] 



Vi 2 \r) 



h 2 m 2 - 1/4 
2M V 2 



The first term in vffl , proportional to m 2 , is the potential which describes the familiar centrifugal force. The negative 
correction term —1/4 gives rise to a centripetal force which from this point on we shall call a quantum anti-centrifugal 
force to emphasize that its binding power arises from quantum mechanics. 

The effect of this contribution stands out most clearly for a particle with zero angular momentum, that is m = 0. 
In this case the effective potential 

y Q (r)^\r)=-^± (5) 
becomes attractive and the radial wave function reads 

ua(r) = VrJa (kr) 

Here and throughout the paper the subscript on the wave functions indicates vanishing angular momentum. 

Hence, the bunching effect of the zeros of Jo takes on a new meaning. It reflects the fact that there is an attractive 
force towards the center x = y = 0. Indeed, the force corresponding to Vq causes a ring-shaped wave packet to spread 
in asymmetric way: It initially spreads faster towards the center than towards the outside [B. 

In three dimensions no such quantum potential exists. We recall || that in this case the effective potential is 
proportional to the quantum mechanical square 1(1 + 1) of the angular momentum quantum number /. Hence, for 
I = the potential vanishes and no asymmetric spreading takes place. 



III. PROPAGATORS 



This dimensionality-dependent spreading is a consequence of the propagators being different in different dimensions. 
In this section we compare and contrast the radial propagators in two and three dimensions. Throughout this section 
we confine ourselves to radially symmetric wave functions, that is, wave functions of vanishing angular momentum. 
This symmetry is preserved under time evolution, therefore it suffices to consider the radial propagator. 

We show that the radial propagator of a free quantum mechanical particle in three dimensions is the difference 
between two propagators of the free one-dimensional radial motion. However, in two dimensions the propagator 
is different, reflecting the attractive quantum anti-centrifugal force. The mathematical reason for this difference is 
exactly the same as in the d'Alembert propagators: It is the difference between a spherical and an ordinary Bessel 
function. 
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A. One dimension 



In this section we lay the foundations for the following sections by briefly recalling the Green's function G^ of the 
free Schrodinger particle of mass M. The Green's function 

G ( g ] (a,* | 5,4 = 0) = Af(t) e M(t)(l - i)! (6) 

with the normalization constant 



and 



allows us to propagate the initial wave function ^Q~\x,t = 0) = $q ( x ) to a later time, such that 

<Z>W(x,t)= dxG^ (x,t\x,t = 0)$ i o ) (x) (9) 



For the remainder of the article it is useful to recall that for t — > the parameter a and the normalization constant 
Af approach infinity. 

B. Two dimensions 

We now use the expression @ for the propagator Gg in one dimension to derive the Green's function of the radial 
motion in two dimensions. For this purpose we take the product of the one-dimensional Green's functions of the 
motions along the x- and the j/-axes together with the initial wave function and integrate over the angular and radial 
variables. We then express the wave functions by the radial wave functions. 

1. Radial propagator 

We start from the generalization 

/OO POO 
dx / dy Gf (x,t\x,0) 4 X) (y,t\y,0) $ 2) ( y/^Tj?) (10) 
- oo J — OO 



of the propagation equation (|J) to two dimensions. Here we have considered a rotationally symmetric initial wave 
function (x, y; t = 0) = "Jq 2 ' (p) , which depends on the radial variable p = \J x 1 + y 2 only. 

When we substitute the expression (||) for the one-dimensional Green's function Gg into ([l(]) and introduce polar 
coordinates p and ip we arrive at 

^\x,y;t) = ¥ ( ?\r;t) 

dpAfe ta(r2+p ^J\f^p- d^e~ aarpcosip ^<i> Q 2) (p) (11) 



The wave function ^ (x, y\ t) depends on the coordinates x and y only through the radial coordinate r = \J x 2 + y 2 . 
Hence, the axial symmetry of the initial wave function is preserved under time evolution. 

We recall the Sommerfeld representation (||) of the ordinary Bessel function Jo, and the propagation equation 
reduces to 

r°° 1 

* 2) (r;4) = / d P Are ia ^ + "^-=y/WSFp2J [2ar P ] Vp^o\p) (12) 
Jo V* 

where we have made use of the definition (Q) of Af(t). With the help of the radial wave functions 
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and 



we find 



t4 2 '(r;t) = VF^ a) (r;<) 
w f ( P )^(p ; t = o) = v^ 2) (p) 



Z" 00 1 
4 2) (r;t) = / dpAAe M(r2+p2) ^=V¥H?p2 J [2arp]4 2) (p) 
Jo V« 



= / dp5 w (r,*|p)^ 2) (p) 
Jo 



This relation suggests the formula 

gP>(r,t | p) = Af(t) e ia ^ r ' + p'^y/Tra{t)rp 2J [2a{t)rp] 



(13) 



(14) 



for the Green's function of the radial motion of a free particle in two dimensions. The key feature of this result is the 
emergence of the ordinary Bessel function Jq. 



2. Alternative form 



It is instructive to represent ||10|] the Bessel function 

MO = I 



in terms of the two Hankel functions 



and 



where we have introduced the abbreviation 

1 



4 e-^-^n^) 



dr- 



(15) 



Wo VCVi-</(20 V^Jo " ^1-^/(20 

We substitute these expressions into the formula ( |l4| ) for the propagator and complete the squares. We arrive at 

g (2) = N e Mr-p) 2 n (2 a rp) - iAfe^+^H* (2arp) (16) 
where we have used the relation 



-in/2 



1/2 



-i-rr/4 



When we compare the radial Green's function in two dimensions (16) to the Green's function Gg in one 



dimension 



we find the representation 

S (2) M I P) = G { s 1] (r,t\ p)H[2a(t)rp]-iG s 1 \r,t\ -p)H*[2a(t)rp] (17) 

Hence the Green's function of the radial motion is the difference of two one-dimensional Green's functions: 
One corresponds to the motion starting at p, the other at the mirror image —p. The two functions reflect the fact 
that the radial coordinate is only defined for positive values. In other words, there is an infinitely steep and infinitely 
high potential well at p — and the radial wave function has to vanish there. This fact guarantees that the radial 
momentum operator is Hermitian for r > 0. 

Moreover, there is a phase difference of it/2 between the two contributions in (p"7|). This phase difference is just 
another manifestation of the interference of plane waves in two dimensions giving rise to an ordinary Bessel function 



as discussed in section II A 



However, the most important feature of ( |17| ) is the fact that each contributing term is the product of the Green's 
function of a free motion and the function Ti. Hence the propagator of a free particle in two dimensions and of 
vanishing angular momentum is not simply that of free motion, but is modified by Ti. 
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3. Short time limit 



We can gain deeper insight into the correction factor TL by considering the short time limit, that is t — > 0. From 
section [n] we recall that in this case a — > oo. The argument TL is the product of a and the radial coordinates r and 
p. Hence we are dealing with a rather subtle limit. 

For r =/= and p^Owe can consider the asymptotic limit of TL for large argument £. When we expand the second 
square root in the definition (051) of TL we find 



2 



OC 



H(0 - — / dre- r (1+*^ 



With the help of the integral relations 

2 

and 



2 /0 Ve-V = i 
o 2 



we arrive at 

We recall the definition (g) of a and the correction factor TL to the free propagator reads 

r ^ n /M \ 1 ?i 2 1 

W*)rp]=«(^pJ =! + ---* 

The second term here is proportional to the quantum anti-centrifugal potential Vq defined in (^|) and the function TL 
takes the compact form 



TL[2a(t)rp] = 1 - l -V Q (^/rp)t ~ exp 



(18) 



This formula shows clearly that for short times the radial motion in two dimensions is determined by the quantum 
anti-centrifugal potential Vq. 

We emphasize that the expression (|l8] ) is only valid for short times and for radial coordinates r ^ and p ^ 0. For 
t 7^ and r = p = the radial Green's function g( 2 ) vanishes, that is 

g^(r = 0,t\p)=g^(r,t\p = 0)=0 

(2) (2) 

as required by the boundary condition u (r = 0) = and Vq (p = 0) = 0. Indeed, this property follows from the 
exact expression (|l4| ) for </ 2 ) in terms of the Bessel function. 



C. Three dimensions 

We now turn to the case of a free particle in three dimensions. Again we consider a wave function \&q (x, y,z;t — 

o) = $< 3) (yi 2 + y 2 + ^ 2 ) 

that at t — is rotationally symmetric, corresponding to vanishing angular momentum, 
that is, I = 0. Time evolution preserves this symmetry. In order to find the Green's function in three dimensions we 
take the product of the Green's functions Gg of the free motions along the x— , y—, and z-axes together with the 
initial wave function, and integrate over spherical coordinates p, and ip. 
We start from the expression 

/>oo />oo />oo 

(x,y,z;t) = dx dy dz 
Jo Jo Jo 

G« (x, t | x)G^ (y, t | y)Gf (*, * I ^ (V^+¥+^) (19) 



G 



for the propagation of the wave function $q 3 ' ) . With the help of the explicit formula (0) for the Green's function 
and spherical coordinates we find 

9®(x,y,z-,t) = 9$\r-,t) = 

•OO p7T fit 

dpJVe la( - r2+p2 *>Af 2 p ddsm8e- l2arpcosf> dipp<S>f\p) (20) 



Since the right hand side depends on the coordinates x, y, and z only via the radial coordinate r = \J x 2 + y 2 + z 2 , the 

(3) 

wave function ^ Q at time t only depends on r. Hence, the radial symmetry of the initial wave function is preserved 
during time evolution. 

The integrand is independent of w and we can immediately perform this and the integration. Indeed, with the 
help of the integral representation (Fy) of the spherical Bessel function jo we arrive at 

/>oo 

^\r;t)= / dpNe %a ^ +p2 HnN 2 pj {2arp)p^\p) 
Jo 

When we introduce the radial wave functions 

u ( 3) (r,t) = r¥ 3 \r,t) 

and 

vi?\p)mug\p,t = 0) = p^(p) 
in three dimensions corresponding to vanishing angular momentum we find the propagation equation 

/•OO 

4 3) M)= / dpNe la{ - r2+p2 HTTN 2 rp ]0 {2arp)vf ) {p) 
dpgW{r,t\p)v<g\p) 

Jo 

This approach yields the Green's function 

g^{r,t | p)=M(t)e la ^ r2+ ^-2a(t)rpj [2a(t)rp] (21) 

i 

for the radial motion in three dimensions. Here we have used the definition (0) of the normalization constant TV. 

When we compare the expression ( |2l| ) for the radial propagator in three dimensions to the corresponding Green's 
function (14) in two dimensions we find, apart from some factors, two changes: (i) In the transition from two to three 
space dimensions we have replaced the ordinary Bessel function J by a spherical Bessel function jo, (ii) moreover, 
the square root of i in the denominator in equation (|l4| ) is now replaced by i. 

These two at first sight minor substitutions have dramatic consequences. As discussed in equation (^) the spherical 
Bessel function is a superposition of an incoming and an outgoing spherical wave with a phase difference of n. When 
we use this representation of jo, we find the Green's function 



g^(r,t \p)=Af(t) 



e ia(t)(r-p)- _ e ia(t)(r+ P y 



for the radial motion in three dimensions. We emphasize that this result is valid for all times. 

Similar to the Green's function (IT^) in two dimensions, the corresponding expression in three dimensions is again 
the difference 

gW(r,t\p)=G ( i\r,t\p)-G < £\r,t\ -p) 

between two propagators of one-dimensional motion: One corresponds to the motion starting at p and the other one at 
—p. However, in contrast to the Green's function in two dimensions, now we deal with a purely free motion represented 
by the Green's function G^: The modification factor Ti of two dimensions is absent since in three dimensions there 
is no quantum anti-centrifugal potential. 

Moreover, there is a phase shift of tt between the two contributions. Its origin is clearly the interference of plane 
waves in three dimensions giving rise to a spherical Bessel function. 
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IV. CONCLUSIONS 



We conclude by emphasizing that the radial propagator of a free particle of vanishing angular momentum moving 
in three dimensions is determined by the spherical Bessel function. It originates from the interference of plane waves 
in three dimensions. Its explicit representation in spherical waves shows that the radial propagator is the interference 
of two propagators of free motion in one dimension. The phase difference between the two contributions is w. 

In two dimensions the radial propagator is determined by the ordinary Bessel function. It originates from interference 
of plane waves in two dimensions. Since an ordinary Bessel function is not a superposition of two cylindrical waves but 
involves correction terms, the propagator is not a superposition of propagators corresponding to free, one-dimensional 
motion. This feature reflects the presence of the quantum anti-centrifugal force in two dimensions. 

In mathematical terms, the dramatic difference between the radial propagators in two and in three dimensions 
results from the different area or volume elements. A comparison between the corresponding equations (|ll|) and j20| ) 
brings this feature out most clearly: In two as well as in three dimensions the phases of the interfering waves, that is 
the phases of the exponential in the angular integration are identical. However, the integration measures arc different: 
In two dimensions we have dip whereas in three dimensions we have d9sm8 = —d(cos9). This difference leads to an 
ordinary Bessel function Jo in two dimensions but to the spherical Bessel function in three dimensions. 

It is interesting to note that this very same reason also causes the difference between the Green's functions of the 
d'Alembert wave equation in two and three dimensions, and the breakdown of Huygens' principle in two dimensions. 
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APPENDIX. GREEN'S FUNCTION FOR THE D'ALEMBERT WAVE EQUATION 

In this appendix we briefly derive the Green's function G^ N ' of the TV-dimensional d'Alembert wave equation 

(22) 



Here r and denote the TV-dimensional position vector and Laplacian, respectively. 

In order to facilitate the comparison between the Green's function of the d'Alembert wave equation and the 
Schrodinger equation of a free particle in the main body of the paper we focus on the cases of two and three 
dimensions. We show that in three dimensions the Green's function of the d'Alembert equation is a delta function. 
In contrast, in two dimensions the Green's function carries a tail behind a sharp edge. We trace this difference back 
to the difference between the spherical and the ordinary Bessel functions. 

Fourier solution in N dimensions 

In order to solve the inhomogeneous wave equation ( [22] ) we make the Fourier ansatz 

G {N) {t,r) = fdk [ d N kG{k ,k)e likoCt -^ 



where k denotes the TV-dimensional wave vector. 

When we substitute this ansatz into ( [22] ) and use the Fourier representation 

of the product of delta functions, we find the algebraic equation 

1 



G(ko,k) = - 



(27T) 



N+l 
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with the solution 



G(k ,k)- {2ll) N + i kl _p {2n) N + i 2 fc 



k 



where k = \ k | . 

Hence, the Fourier transform G enjoys simple poles at fco = ±/c. Different paths in the complex plane fco define 
different Green's functions. Since in the present paper we are only interested in the dependence of the Green's function 
on the dimensions we choose one particular propagator, namely the retarded propagator G^et . 

The retarded Green's function is defined by a path just below the real fco-axis. For KOwe can close the path by 
a circle at infinity in the lower complex plane without changing the value of the integral. Since in this case we do not 
include any poles we find 

G^(t,r)=0 

On the other hand when t > we can close the path in the upper complex plane without changing the value of the 
integral. The residue theorem 

dz — - = 27ri/(zo) 
z-zo 

where the integration path circumvents the pole zq in the counter-clockwise direction allows us to perform the 
integration over ko- We find the expression 

(JV) _ (~j) f d N k , ikct ikct\ -ik-r 



^ ct " (2n) N J 2k [C 6 > 

We can combine the results for positive and negative values of t in the form 

«(*,rl = 6(t)^ J d-k^e'^ (23) 

It is instructive to cast this formula into a slightly different form by separating the iV-dimensional volume element 



d N k^k N - i dkd N ~ i 



into the integration over the length of the iV-dimensional wave vector and the N — 1 angles contained in dui. 
The Green's function then reads 



G'Z> (t, r) - Q(t) ^ dkk" - 2 S in(kct) 



toe 



The retarded Green's function of the iV-dimensional d'Alembert wave equation is determined by a two-fold interfer- 
ence: (i) the interference of plane waves in the N — 1 dimensional space of angles and (ii) the interference of these 
expressions due to the integration over all wave numbers. In (i) all plane waves have the same amplitude and the 
same wave number. However, they differ in their propagation direction. All directions appear with equal weight but 
the surface element can put different weight on different directions. It is this contribution that creates the difference 
between the propagator in two and in three dimensions, as we show now. 



Two dimensions 

We first consider the case of two dimensions, that is N — 2 and perform the remaining integrations over the 
two-dimensional wave vector k in polar coordinates. With the help of the area element 

d 2 k = dk k dip 

and the Sommerfeld representation (0) of the Bessel function Jo, the Green's function 
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1 f°° 1 f n 

G^ = e(t)— / dksin(kct) — 

27T JQ 27T J— it 

following from (23) takes the form 



dip e 



— ikr cos if 



The integral relation 



finally yields the expression 



i r°° 

GlH=@(t)— dksm(kct)J (kr) (24) 
2tt Jo 

Q(a 2 -b 2 ) 



dtsm(Za)MZb) 



Va 2 - b 2 



2tt ^/(ct) 2 - r 2 

for the retarded Green's function of the d'Alembert wave equation in two dimensions. 

This Green's function enjoys a sharp wave front due to the Heaviside step function. Moreover, it is followed by a 
decay which obeys an inverse square root dependence. 

Three dimensions 

We now turn to the three-dimensional case where the volume element in wave vector space reads 

d 3 k = dk k 2 d9sin9 dip 
The retarded Green's function then takes the form 

i y* OO />7T -| |>7T 

Gill = e(t)---^ dksm(kct)k d6 sin 6e- lkrcos 9 — dp 
(2tt) 2 J J 2tt J_ 7I 

which after simple angular integrations with the help of the integral representation (Q) of the spherical Bessel function 
jo reduces to 



dk sin(fcci) k jo(kr) 



When we compare this expression to the corresponding formula ( |2~4| ) for two dimensions we recognize that the ordinary 
Bessel function Jo has been replaced in three dimensions by the spherical Bessel functions jo. This has important 
consequences for the remaining integration over k. Indeed, when we use the representation (|l]) of the spherical Bessel 
function as a sine function we find 

GUI = 6(0 / dksin{kct) sin(fcr) 

nr 2n J Q 

which with the help of the trigonometric relation 

sin a sin /3 = ~ [cos(a — (3) — cos(a + /?)] 

yields 

GUI = Q(*)^— 7T / dk{cos[k{ct - r)] - cos[k(ct + r)]} 
Ziir Ztt Jo 

or 

Gill = Q(t)^-[S(ct-r)-S(ct + r)} 

Due to the Q function, the retarded Green's function is only non-vanishing for t > 0. Hence, only the first delta 
function makes a contribution and we arrive at the expression 

Gill - e(t)*^ 

ret 47rr 

for the retarded Green's function of the d'Alembert wave equation in three dimensions. The integration of the spherical 
Bessel function instead of the ordinary Bessel function has created the delta function in the propagator instead of the 
square root and function. 
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